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Abstract —There are now several works on the use of the 
additive inverse Gaussian noise (AIGN) model for the random 
transit time in molecnlar communication (MC) channels. The 
randomness invariably causes inter-symhol interference (ISI) 
in MC, an issne largely ignored or simplified. In this paper 
we derive an upper bound and two lower bounds for MC 
based on amplitude shift keying (ASK) in presence of ISI. 
The Blahut-Arimoto algorithm (BAA) is modified to find the 
input distribution of transmitted symbols to maximize the lower 
bounds. Our results show that over wide parameter values the 
bounds are close. 

Index Terms —Capacity, Inter symbol interference, lower 
bound, molecular communication. Upper Bound, AIGN channels. 

I. Introduction 

Molecular communications (MC) systems encode informa¬ 
tion in the concentration HI, time of release El, the number 
of molecules released in a time-slot and the type and ratio of 
molecules 13 . This technique has potential applications in en¬ 
vironments where electromagnetic waves cannot be used, e.g., 
in buried pipelines or in medical applications with embedded 
devices. This paper considers the case where information is 
encoded in the number of molecules released in each time- 
slot, i.e., molecular amplitude shift keying (ASK). 

The propagation of molecules in a fluid medium is governed 
by Brownian motion a, possibly with a drift velocity El. 
In 121 such propagation was analyzed and it was shown that 
the molecules experience a propagation delay that follows an 
inverse Gaussian distribution. Based on this so-called inverse 
Gaussian noise (IGN) channel, in El, 121 expressions and 
bounds for channel capacity are derived when information 
is encoded on the release time of molecules. The capacity 
of MC when information is encoded on the concentration of 
molecules is studied for binary communications in l6l, |2l and 
for binary and 4-ary communications in |3. In these works, 
the aggregate distribution of the number of arrived molecules 
in a time slot is approximated by a Gaussian distribution. 

In molecular ASK modulation, molecules diffuse to the 
receiver where they are detected and removed from the system. 
One of the main challenges in such a diffusion-based com¬ 
munication system is the resulting inter-symbol interference 
(ISI); due to the random propagation time, molecules may 
arrive over many time-slots. Of interest here, therefore, is a 
capacity analysis of a more likely system which suffers from 


ISI. In a, disregarding the presence of ISI, a binary ASK 
scheme by considering the life expectancy of molecule with 
AIGN model for propagation time is studied (output symbols 
are independent). Indeed, if the system does not suffer from 
ISI, in model wherein molecule release and detections are 
perfect, communications would also be perfect. In (91, the 
capacity and the probability of error of binary and 4-ary ASK 
schemes are investigated in presence of ISI. However, for ca¬ 
pacity analysis, researchers consider maximization of symbol- 
by-symbol mutual information, which is a lower bound for 
li.i.d or the achievable information rate for ISI channel with 
independent, identically distributed (i.i.d) inputs ITOl . 

Capacity analyses of conventional ISI communications has 
a strong presence in the literature. Starting from the seminal 
work in El, several methods have been proposed to derive 
the capacity of an ISI channel. The problem of determining 
the channel capacity is completely solved in the case of an 
unrestricted input distribution but remains open for distribu¬ 
tions over a discrete input alphabet lfT3 . Several lower and 
upper bounds to the achievable rate with discrete inputs and 
Gaussian channels have been obtained in the literature m, 
ca.A simulation-based approach is applicable for calculating 
capacity in the case of i.i.d. inputs with binary modulation 
ifTrll . ifTSll . This approach is based on a trellis whose number 
of states corresponds to the memory of the ISI channel. 

In this work, we consider a general molecular ASK scheme 
in presence of diffusion-induced ISI. Time is slotted and the 
information is encoded in the number of molecules released in 
a time-slot, and the receiver counts the number of molecules 
received within each time-slot. We obtain the probability of 
molecules arriving within a specihc time-slot using the AIGN 
model, in turn leading to a binomial distribution on the number 
of received molecules in each time-slot. We propose two lower 
bounds and an upper bound for the Ii,i.d of such a channel, 
which is, here, called capacity. To provide a tractable analysis 
we restrict the effect of ISI to one time-slot, i.e., we assume 
the molecules that do not arrive within two time-slots have 
disappeared. The hrst lower bound is based on a symbol-by- 
symbol maximization of mutual information, while the second 
lower bound uses information of next received symbol for 
decoding. The lower bounds are maximized by optimizing the 
input distribution, which in turn closes the gap between the 
lower and upper bounds. 


The remainder of this paper is organized as follows, the sys¬ 
tem under consideration is presented in Section II. The channel 
characteristics including the effect of ISI are developed in 
Section III. In Section IV an upper bound and two lower 
bounds of MC with ISI are proposed including an optimization 
of these bounds. Section V presents some numerical results; 
finally Section VI wraps up the paper with some concluding 
remarks. 


II. System Model 


In our system, the transmitter is a point source of identical 
molecules. At the beginning of every time slot, with length of 
T, it conveys a message by releasing X molecules into a fluid 
medium. Here, X with 0 < X < is a random variable 

and Xiaax is the maximum number of molecules released in 
any time slot. The transmitter therefore uses amplitude shift 
keying (ASK), i.e., the message is encoded in the number of 
molecules released. Once released, the transmitter does not 
affect the propagation of the molecules. 

Molecules propagate between the transmitter and receiver 
by Brownian motion characterized by a diffusion constant d 
and (positive) drift velocity v. At the receiver, all received 
molecules are absorbed and removed from the system. Impor¬ 
tantly, for our analysis, we assume that a molecule transmitted 
in time-slot m either arrives in the same time-slot, the next 
time-slot (m-fl) or disappears. Essentially, for tractability, our 
analysis focuses on ISI within one time-slot. We assume that 
everything else within the system operates perfectly, i.e., the 
only randomness in our model is the propagation time (causing 
a randomness in the number of molecules that are received). 
Using the AIGN analysis in El, the cumulative distribution 
function (CDF) of the propagation time is given by 


Fw (w) 



Here $ (■) is the CDF of a standard Gaussian random variable; 
also, ^ = l/v and A = Z^/cr^ in which I is the distance between 
transmitter and receiver and tr^ = d/2 is the variance of the 
associated Weiner process El and d is the diffusion coefficient. 

Let qk denote the probability that a molecule arrives within 
the fcth time-slot following its transmission, i.e., a molecule 
transmitted at time 0 arrives in time interval {{k — 1)T, kT], 
where T is the time-slot duration. We have qk = TV (kT) — 
Fw {{k — 1) T). The probability of a molecule arriving in the 
same time slot that it is released is qi = Fw (T). 

The probability of receiving Vi = y molecules at receiver, 
when transmitter release Xm = x molecules is binomial: 


PTiY^=y\Xm 



y ^ 9i(l - qiT ^,0 < y < T 

0 ,y<0,y>x 

( 2 ) 


We denote by a^;, 1 < x < ATmax, the probability of the 
transmitter releasing X^ = x molecules. While (|2]) ignores 
the effect of ISI, we consider its impact in coming sections. 

HI. Communications in the Presence of ISI 

With ISI we must consider transmit/receive sequences and 
the average mutual information per channel use is given by 

I{X^-Y^+'^)= lim , (3) 

L—>-oo L ^ ' 

where, X^ = [Xi,...,Xl] and = [Fi,...,V+fc] , 

denote input and output sequences of length L and L+k evalu¬ 
ated for a given input distribution Px^ (X^). This determines 
the achievable rate of reliable communication through this 
channel with this specific input distribution and the channel 
capacity is the supremum of this mutual information over all 
allowed input distributions. This is in contrast to the case of 
a discrete memoryless channel (DMC) without ISI and all 
molecules either arrive within the time-slot of transmission 
or disappear. We denote the maximum value of the mutual 
information of this DMC as Cdmc- 

For tractability, we will focus on the case of independent 
and i.i.d. inputs, i.e., Pxl — a x ... x a. = a^, where a = 
[oo,..., uXmax]’ ^ simplification suggested by the work in ifTOll . 
We denote the resulting average mutual information as Ii,i,d. 
Note that this mutual information is a lower bound of the 
expression in ©, but with a slight abuse of notation we will 
call the resulting mutual information the “channel capacity” 
given by 

C= lim \sn^l{X^-Y^+^) . (4) 

L—]_j 

Despite this significant simplification, calculating the chan¬ 
nel capacity in (IHi appears intractable and, in this paper, we 
develop lower and upper bounds on this expression. 

At the start of time slot k, Xk € {0,1,..., Ai„iax} molecules 
are released. We begin by considering the interference from 
k — 1 previous time-slots. Let Pk (n) denote the probability 
that n molecules arrive from the k — 1 previous time-slots, 
in time slot k. Furthermore, let Zk denote the number of 
transmitted molecules (at time-slot k) not received in the 
same slot and let Nk-i denote the total number of inter¬ 
fering molecules from the previous {k — 1) time slots, i.e., 
Nk-i G {0,1,..., (fc - l)Xi„ax}- Finally, letting Yu denote 
the number of molecules received in time-slot k we have 

Yk = Xk — Zk + Nk-i, ( 5 ) 

we have Yk G {0,1,..., Aimax^}. The transition probability 
when Xk = x molecules are transmitted, it is given by 

PY^iXk (.y\x) = 

(1 - qi)''Pk (y) + ^ ^ ^ '?i(l - qiT~^Pk (y - 1) + ... 

( a; % ) (1 - qi) Pk (y - (t - 1)) + qi^Pk (y-x), 

y — 0: (^max . 


(6) 


Probability of receiving y molecules for k — 1 time slots, 
Pk (y), y = 0,1,Xinaxfc- , simply Can be calculated 

by induction. Hence we have 

Pk {y) = (ao + + (1 - Pk-i [v) 

+ (quai + ...+ ) qk{l - qkf---\x^^}j X 

Pk-i (y ~ 1) + 

+ 1 yfe OX,aax-l + I - 1 ) ^ 

®Xmax) Pk -1 (y ~ (^max ~ 1)) + 
q^ ““OjCmax^fc-l iy ~ ^max) ) 

(7) 

where = Pr (xi = i) for i G {0,..., Xmax}- Using induc¬ 
tion we can compute Pk (y). In general, we can write 

Pk (0) = (ao + ... + (1 - qkf-'ax^^}j Pk-i (0), ( 8 ) 

Pk (1) = (ao + ... + (1 - qkf-'ax^^}j Pk-i (1) + 

J ^gfcai + ...+ (^ 
a-x^^JPk-i (0), 

(9) 


Pfc((fc-1)-Ymax) = 

(ao + ... -|- (1 — qk) ”“ax„ax) Pk-i {{k — 1) ^max) + 

J ) Qkai + ...+ ^ qk{l - X 

^-’fmax) Pk-1 {{k — 1 ) Xmax ~ 1) + 

®^max) Pk-1 {{k — 1 ) X„iax ~ 2 ) -f 


qk) X 

( 10 ) 

Pfc ((fc - 1) X„,ax + 1) = 0, (11) 

IV. Upper and Lower Bounds 

We now provide two lower and an upper bounds for MC in 
presence of ISI. While the expressions in (| 6 ]) and O are valid 
for any k, for tractability, we now assume that the ISI only 
affects the next time-slot, i.e., molecules are received within 
two time slots or disappear. This is equivalent to fc = 2 and 
we have Pi( 0 ) = 1 , Pi (n 7 ^ 0 ) = 0 . 

A. Lower bound 1 

In the first lower bound we consider the effect of ISI 
degradation on mutual information between input and output 
symbols. This is a lower bound on the channel capacity 
because this measure ignores the memory; essentially, we 


^Wmax 1 I 

% OW^ax-l + 


Xmax 
X^nia.x 1 


.--^max 1 


% 


( 1 - 


OAT^ax) Pk-1 {{k — 1 ) Xmax — (-^max — 1)) + 




^X^axPk—1 {{k 1 ) Xmax X^ax) 5 


consider a DMC but with an additional source of measurement 
error due to molecules from the previous timeslot. Indeed, 
comparing the capacity of the DMC (Cdmc) this lower 
bound measures the error of ignoring ISI. This lower bound 
relates to lower bound of h.i.d in discrete Gaussian channel 
with ISI in Eo). Hence we have 


Ilb, = I {Xm-, Urn) = H (Y^) - H{Y^\X^) 

2-^max 

= - E P(ym)log(p(ym)) + 

ym —0 

- ^ma x a;m-|-A:,„ax 

E E p(ym|a;m)logp( ym I ^m) ; 

Xrn—0 ym —0 

where p{ym\xm) is given by 


( 12 ) 


p(yn 


{i-qir--\\x 


0 = < 


E 

i=0 

^max 

E 

(1 - <Xm+X^ 

^^ym P Xm -f Xn 


J 

ym - i 

j-iVm-i) 


Vm —1 


(13) 


By averaging over a;™ on p (y^l a;^) , P (ym) is given by 


- ^ma x ^771 

P{Vm) = E E 

Xm—0 i—0 

-^max / q 

S iv -i 

j^y-i \ ym ^ 


(l-yi)"ylx 


ajqr~\^ - 92> 


i-iVm-i) 


(14) 


B. Lower Bound 2 

We assume that molecules are received within two time- 
slots or disappear; hence the transmitted symbol in current 
time-slot only affects received molecules in current and next 
time-slot. We consider mutual information between transmit¬ 
ted symbol in current time-slot and received symbols in current 
and next time-slot 

kLB2 — 7 (Xm—17 Ym— 1, Ym ) 

= H {Ym,Ym-l) - H {Ym, Ym-l\ Xm-l) ,, 

(b) 

= P (Xm- 1 ) + 77 ( Xm I Xm- 1 ) - 77 ( Xm- 1 1 Xm- 1 ) 

-7T (Xml Xm-l, Xm-l) 7 


where (a) is obtained based on definition of mutual information 
and (b) based on definition of joint entropy. We consider the 
channel in steady state regime, hence P (ym-i| a^m-i) = 
P ( Pm I Xm'), which is given in (O. Also, p (ym| ym-i, a^m-i) 
is given by 

p{ym\ym-l,Xm-l) = 

E P ( I ym— 1 7 l) P ( ym I ym— 1 7 ^m —l) MAt 

(b) Amax 

E P (^m) P ( ym I Pm -17 ^m, ^m— 1 ) 7 

Xm—0 


where (a) is obtained based on law of total probability, (b) is 
obtained based on independence of Xm from Pm-i and ccm-i 
which is due to the i.i.d assumption of the input distribution 




and causality. By averaging over Xm-i in (fT^ . p{ym\ Um-i) 
is given by 


P{y7n\Vm-l) = 

- ^ma x 

X; P{x^_i\y^_i)p{ym\yni-l,Xm-l) 

-^max -^max 

= Yj p{Xm-l\ym-l) Y p{xm)^ 

a^m-l=0 Xm—0 

v(v \v ^ X X 1 i p(ym-lkm-l)p(a^m-l) 

P\ym\ym—l'i ^rrii'^m — 1) 2^ p(y^_i) ^ 

- ^ma x 

P {p^m) P ^ym\ym —It‘^rri’^• 

Xm—Q 


(17) 


where (a) and (b) are obtained based on law of 
total probability, and Bayes’ rule, respectively. Also, 
P (ym I ym —15 Xrm Xm—i) is giveu by 


P ( Pm I ym — li Xmi 1) 

(' 




P ( y m—2 I —l) ^ 

y'm-2=0 

P ( ym I ym —15 Xjn^ Xm—l ; y m — 2 ) 

(b) / / I t 

= L p(.y m- 2 | ym-l,a:m_ijx 

y'm.-2=0 

P ( ym I ym —15 ^mj ^m—1; y m— 2 ) 

p(yTn-i,a;m-i) 

P ( Pm I Pm — li Xrm Xm—l ; P m — 2 ) 

p(ym-l|p',n_2,3^m-l)p(a:m-l|y'm-2)p(y'm-2) , 

p(yTn-l|a;m_i)p(a;m_i) 


-^max 

= E 

^'^-2=0 


(£) 

y'm.- 2=0 

P (pml Pm-1) 3^m—Ij P m— 2 ) 

p(t/m-l|y'„i_2,a:m-l)p(a:m-l)p(y',n_2) , 

p(pm-l|a;m-l)p(a:m,-l) 


(e) 


= E 

y'm-2=0 

P ( Pm I Pm-1; ^m j Xm— 1; P m— 2 ) 

p(Pm.-l|p'm-2.^m-l)p(y',n-2) , 

p(pm-i|a;m-i) 


= E 

P'm-2=0 

P ( Pm I Pm — li Xmi Xm—l iP m — 2 ) 

(_£) p[x'm-l\xm,-l)p{y'rn- 2 ) ^ 

y m- 2=0 

P ( Pm I Pm —1 j 3^m—1 j P m — 2 ) ; 


(18) 


where (a) is obtained based on law of total probability and 
Pm -2 i® number of received molecules at the end of time- 
slot m — 1 from transmitted molecules in time-slot m — 2; 

(b) is obtained due to the independence of p ^_2 from Xm', 

(c) is obtained based on Bayes’ rule; (d) is obtained based on 
the joint probability formula, (e) is obtained based on indepen¬ 
dence of ym -2 Xm-i and the joint probability formula; (f) 
is obtained due to the fact that Pm-i = x'm-i + ym- 2 ’ where 
x'm-i denotes the number of absorbed molecules at time-slot 
m — 1 at end of time-slot m — 1. Based on dehnition of y ^_2 


we have 


Amax / T ri 

p(y'm- 2 )= E «.,.-2 

Xm-2—0 \ y m — 2 

(l_q2) — 


(19) 


P ( ym I ym —It ^m i ^m—1 1 y m— 2 ) 

P yPml ^mt ^ m — l~ ^m — 1 (ym—1 P 771 — 2 )) “ 

{ ■^ ( ^m ^ xaim —* i f ^ m — 1 

.=..E-A 

qY *( 1 - 92 )'^"*"' '‘\ym < Xm + x”m -1 

H ’Pm ^ Xm “b X m—h 

( 20 ) 

(a) is obtained because tCm-i = x'm-i+x'm-i, where x'm-i is 
the number of remaining molecules at end of time-slot m—1 
from transmitted molecules in the same time-slot. Moreover, 
based on dehnition of x'^_i , P (a:'m-i| Xm-i) is given by 



P {x I Xm — l') — 


Xm—l 
X m—1 


9i 


‘(1-Pi) 


Xm — l—X m — 1 


( 21 ) 


C. Upper Bound 

We assume that molecules are received or disappear after 
two time-slots. Hence, if we send symbols and wait two 
time-slots before transmitting the next symbol we have an 
interference-free channel. This is equivalent to the DMC case 
with the binomial transition probabilities of (l2]i where qi is 
replaced by qu — Fw{2T). The mutual information of this 
channel is concave which can be maximized using the Blahut- 
Arimuto algorithm (BAA) ESI. We therefore have an upper 
bound as follows 


IjjB — mux I (2fm ’ l^m) — mux [77 (Ym) H ( Ym I -Am)] 1 

a a 

( 22 ) 


where P{ym\xm) is given by 


P ( Pm I Xm) — 


/ Xm \ ym \Xm—Vm ^ 

„ % (1 - <lu) iVm < Xn 

\ ym J 


iPm P Xm- 


(23) 


and qu — Fw (2T). This upper bound relates to ’’matched 
hlter bound” of discrete Gaussian channel with ISI in IfTOl . 


D. Matched Filter 


We now provide the achievable rate when using a matched 
hlter. We transmit symbols and wait two time-slots before 
transmitting another symbol. We again have an interference- 
free channel. Naturally, although this scheme may be appeal¬ 
ing in a bit per transmission (channel use) sense, it pays a toll 
when we consider the rate (in bits per second) as it uses two 
time slots. By this dehnition. Imp is given by 


Imp = I {Xm-liYm-l,Ym) = H (Ym, Ym-l) - 


H{Ymi Tm-11 Xm-l) = H (Ym-l) + H (Ym\ Ym-l) " 
H(Ym-l\Xm-l) - H(Ym\Xm-l, Ym-l) , 


(24) 








where (a) is obtained from the definition of mutual information 
and (b) from the joint entropy formula. We measure the 
number of molecules at end of each time-slot. In this case 
P {ym-i\xm-i) is given by 

P ( ym—l I ^m —l) — 


ym — l 


0 


ym—l ^ ^m —1; 

j^m—1 ^ ^m—1- 


Also, p{ym\ym-i,Xm-i) for ym -1 < Xm -1 is givcii by 

P ( ym I —li ^m—1) — 

Xm-l - ym-1 \ ^ 

Vm J 


(25) 


(26) 


Pm ^ Xjji — \ ym—l 

ym ^ ^m —1 ym — l) 


0 


Using law of total probability on (l26] t we have 

-^max 


p(ym|ym-l)= I] p{ym\ym-l,Xm-l)p{Xm-l\ym-l) 

^m-l=0 

(^) , I . p(y^_i|x^_i)p(x^_i) 

/ V P \ ym\ym—1^ ^m — 1 J p(ym-i) ’ 


^m-l=0 


(27) 

where (a) is obtained based on Bayes’ rule, also. Imp is 
calculated by derived a from optimizing Cjjmc- Clearly, 

(b) ( 2 o) 

I {Xm-l',Ym ) > I (Xm-l, Ym-l) = IpBi, 

where (a) is obtained from chain rule in mutual information 
and (b) is obtained based on non-negativity assumption of 
mutual information. Hence, we have following result 


IlBi < IlB2 < li.i.d < luB- 


E. Optimizing The Lower and Upper Bounds 

We are not able to show the concavity of the lower bounds 
IlBi and Ilb 2 with respect to Qx^', however, we can modify 
the BAA ifThl to find a local maximum, and as a result, we can 
calculate these two lower bounds. Each element of channel 
transition probablity matrices pU^h)^ ^ g {1,2} for lower 
bounds 1 and 2 are given by 

=[PiYm=ym\Xm=Xm)], (29) 

[p (Ym Vm,Ym—l ym—l|2fm—1 Xm — 1 

(30) 

The size of and PUB 2 ) (Xmax + 1) x 

(2Xmax + 1) and (Xmax + 1) X (2Xi„ax + 1)^ respectively. 
Also, ym,ym —1 C [0,..., 2Xjnax] and Xm € [0, ...,Xaiax] ■ For 
matrix Q with size of pUBn)^ g 2 } , let 

j(a,p(^^'*),Q) =y] y] a^pj^^'hog^. (31) 

V / 3 ^ dj 

where denotes the element (ijj)'^ of Q. Then the follow¬ 
ing is true 

1) /lb,, = maxmax J (a, Q) . 

a Q ^ ^ 


2) For fixed a, J (a, Q) is locally maximized by 

„ -pUBn) 

Qlj = ' , ■ (32) 

3) For fixed Q, J (a, Q) is locally maximized by 


exp(EzP5^^'‘^logQi,i) 

=-S-FTTTl-^ (33) 

E.exp (E.Pjf'-^iogQL.) 

The algorithm iterates between Oj derived in ( 1331 ) and the 
transition probability matrices pUBh) jjj ( |2^ or (|30] |. This 
procedure is repeated until the convergence of aj. 

In contrast to the lower bounds. Tub is. a concave function 
in terms of Oj. Hence using the standard BAA ifT^ the upper 
bound can be maximized. 


V. Numerical Results 

In this Section the bounds of information rate of ASK 
modulation in channel with ISI and AIGN model for trans¬ 
mission time are evaluated and compared with each other. 
Imp and the Cdmc- Note that in the DMC case, molecules 
arrive within their time-slot or disappear; as a result the 
corresponding performance curves are simply provided for 
insight and context on those of the ISI channel, as opposed to 
a direct comparison. We also study the effect of parameters in 
molecular medium such as I, v and 

Fig. 1 plots Ilb 2 and IpBi as a function of T with v = 1, 
I = 10“^ and cr^ = 1 for the optimized input distribution 
from (l33l l and a uniform distribution. It can seen that using the 
optimized distribution, the bounds Ilb 2 and IpBi is improved 
significantly. Moreover using the optimized distribution of (l33t 
increases the difference between Ilb2 and IpBi to increase 
in comparison with the uniform input. 

Fig. 2 plots Ilb2, IlBi, Imp, Cdmc and Iub versus 
T for different values of I with u = 1, cr^ = 1 and 
Xmax = 7- It can be observed that by decreasing I all bounds 
increase and converge to log 2 (X^ax + 1) = 3 , which is 
the entropy of the source. Crucially, over wide ranges the 
upper and lower bounds are close. Also, due to reduced effect 
of ISI, by reducing I the gap between all derived bounds is 
reduced. Furthermore, Since by increasing T the probability 
of receiving molecules in current time slot is increased. Imp 
also converges to Cdmc for larger values of T . 

Fig. 3 plots Ilb 2 , IlBi, Imp, Cdmc and Iub versus T for 
different values of cr^ with v = 1 , I = 10“^ and X^ax = 7- 
All capacity bounds are increasing functions of a. While this 
seems counter intuitive, an increasing a causes reduced ISI - 
this is consistent with this is consistent with the fact that qi 
is an increasing function of a for low-to-medium values of v, 
i.e., the increased variance in position increases the probability 
of the molecules arriving with two time-slots after release. For 
large u, ~ 1 and ISI is not really an issue. It is worth noting 
that when information encoded in time-of-release as in ||2l, the 
mutual information is not monotonic in a. 







Fig. 1. IlB 2 ^LSi ill terms of T with i; = 1, / = 10 = 1, 

Xmax = 7 for optimized input distribution and uniform distribution. 



Fig. 2. IlB 2 ^ ^LBi, Cdmc and luB in terms of T for different 

values of I (black curves, I = 10~^, blue curves, I = 10~^, red curves, 
I = 10“"^) with V = 1, = 1, Xmax = 7. 


Fig. 4 plots Ilb2^ IlBi. Imf, Cdmc and Iub versus 
T for different values of v with — 1 ^ I = 10 “^ and 
-^inax — 7 . As in El, increasing drift velocity increases 
mutual information due to reduced ISI. 

Comparison the results of the three recent hgures, we note 
that the bounds are most sensitive to the transmitter-receiver 
distance I and drift velocity v while not being as sensitive 
to the diffusion constant tr. Our numerical and simulation 
studies based on performance of an ML (maximum likelihood) 
detector in different settings (not reported here), suggests that 
indeed the presented ISI model with one time-slot memory is 
viable over a wide range of time slot durations, T. 

VI. Conclusions 

In this paper we consider ASK-based MC with ISI. Specif¬ 
ically, we derived two lower bounds and an upper bound on 
capacity (albeit with the signihcant simplification of i.i.d. in¬ 
puts). Our results showed the lower bounds are improved using 
an optimized distribution for the input probability of symbols. 



Fig. 3. IlB2’ IlBi, Imp, C'omc and Iub in terms of T for different 
values of cr^ (black curves, cr^ = 1, blue curves, cr^ = 10, red curves, 
(T^ = 100) with tl = 1, Z = 10“^, Xmax = 7. 



Fig. 4. IlB2^ IlBi, Imp, Cdmc and Iub in terms of T for different 
values of v (black curves, u = 1, blue curves, v = 100) with I = 10“^, 
= 1 and Xmax — 7. 

Also, our results showed that over wide parameter values the 
lower and upper bounds are close. The analysis presented here 
is useful in understanding the sensitivity of a MC system to 
the various parameters. 
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